Light quark energy loss in strongly-coupled AT = 4 supersymmetric Yang-Mills plasma 
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We compute the penetration depth of a light quark moving through a large A^c, strongly coupled 
A/' = 4 supersymmetric Yang-Mills plasma using gauge/gravity duality and a combination of analytic 
and numerical techniques. We find that the maximum distance a quark with energy E can travel 
through a plasma is given by AXnaax(^) = (C/T) {E/T^/Xf^'' with C ^ 0.5. 



I. INTRODUCTION 

The discovery that the quark-gluon plasma produced 
at RHIC behaves as a nearly ideal fluid [1] [2] has 
prompted much interest into the dynamics of strongly 
coupled plasmas. Hard partons produced in the early 
stages of heavy ion collisions can traverse the resulting 
fireball and deposit their energy and momentum into the 
medium. Analysis of particle correlations in produced 
jets can provide useful information about the dynam- 
ics of the plasma including the rates of energy loss and 
momentum broadening [3l H], as well as the speed and 
attenuation length of sound waves [4]. 

Gauge/gravity duality El [71 [8] is a useful tool for 
the study of dynamics of strongly coupled non-Abelian 
plasmas. Although no gravitational dual to QCD in 
known, gauge/gravity duality has provided much insight 
into the dynamics of various theories which share many 
qualitative properties with QCD. The most widely stud- 
ied example is that of strongly coupled A/" = 4 supersym- 
metric Yang-Mills theory (SYM). The deconfined plasma 
phases of QCD and SYM share many properties. For ex- 
ample, both theories describe non-Abelian plasmas with 
Debye screening, finite spatial correlation lengths, and 
long distance dynamics described by neutral fluid hy- 
drodynamics. When both theories are weakly coupled, 
appropriate comparisons of a variety of observables show 
rather good agreement [9l [TOl [11]. This success, com- 
bined with the lack of alternative techniques for study- 
ing dynamical properties of QCD at temperatures where 
the plasma is strongly coupled, has motivated much 
interest in using strongly coupled A/" = 4 SYM plasma 
as a model for QCD plasma at temperatures of a few 
times Aqcd (or 1-5 Tc ^ T < 4Tc). (See, for example, 
Refs. [T2l|l3|ll[T5l[lIl[IZl[li[Tl[ia 
and references therein.) At least for some quantities, this 
is quite successful. In particular, the value of the shear 
viscosity to entropy density ratio [El [26], r]/s = l/47r, in 
strongly coupled SYM is in rather good agreement with 
estimates which emerge from hydrodynamic modeling of 
heavy ion collisions [27j. 
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In the limit of large Nc and large 't Hooft coupling 
A = g'^Nc, the gravitational dual to AT = 4 SYM is de- 
scribed by classical supergravity on the ten dimensional 
AdS^ X geometry [5^. Studying the theory at finite 
temperature corresponds to adding a black hole to the ge- 
ometry |6]. The corresponding A dS-S chwarzschild (AdS- 
BH) metric is given below in Eq. (2.1 ). Fundamental rep- 
resentation quarks added to the JV = 4 theory theory are 
dual to open strings moving in the lOd geometry. In the 
limit of large A, where the string action and the energy 
both scale like a/A, quantum fluctuations in the string 
worldsheet are suppressed and the dynamics of strings 
may be described by the classical equations of motion 
which follow from the Nambu-Goto action. 

The dynamics of strings corresponding to heavy quarks 
have been intensely studied by many authors. The en- 
ergy loss rate for heavy quarks moving through a SYM 
plasma has been studied in Refs |13l [141 128J |29J , and the 
wake produced by a moving heavy quark was computed 
in Refs. [30l [SH [321 [33] . 

Analogous studies for light quarks have yet to be com- 
pleted. In Ref. [34 the charge density of massless quarks 
moving through a SYM plasma was studied, and it was 
shown that there are string states which are dual to long- 
lived excitations (ie., quasi-particles) in the field the- 
ory. In particular, the charge density of highly energetic 
light quarks can remain localized for an arbitrarily long 
time, and can propagate arbitrarily far before spread- 
ing out and thermalizing. In Ref. an attempt was 
made to estimate the penetration depth of a gluon mov- 
ing through a strongly coupled SYM plasma. The results 
of Ref. [3F were obtained by assuming that the endpoint 
of a (folded) string follows a light-like geodesic in the 
AdS-BH geometry; full solutions to the string equations 
of motion were not obtained. The authors of Ref. [35] 
tried to roughly characterize the relationship between the 
string's energy and momentum and the parameters of 
the geodesic, and suggested that the maximum distance 
a gluon of energy E can go before thermalizing should 
scale as Axmax E^^^. The same scaling relation has 
also been discussed for i?-current jets in Ref. [36 . 

Although the estimates made in Ref. [35] are generally 
plausible, we believe that it is clearly desirable to per- 
form a quantitative, controlled study of the penetration 
distance of light quarks (or gluons) in a strongly coupled 
plasma. This is a key aim of this paper. 

It should be emphasized that we are concerned with 
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studying the propagation through the plasma of ener- 
getic excitations which resemble well-collimated quark 
jets. The open string configurations we consider may 
be regarded as providing a dual description of dressed 
quarks, with high energy, moving through a non-Abelian 
plasma. We are not studying the result of a local current 
operator acting directly on the strongly coupled A/" = 4 
SYM plasma. (See, however, Ref. [37 .) Our motiva- 
tion is similar to that of Ref. [17], in which weak cou- 
pling physics in asymptotically free QCD is envisioned 
as producing a high energy excitation, whose propaga- 
tion through the plasma is then modeled by studying 
the behavior of the same type of excitation in a strongly 
coupled A/" = 4 SYM plasma. 

The energy loss rate for a heavy quark depends only 
on the quark's velocity, the value of the 't Hooft coupling 
A, and the temperature of the plasma through which the 
quark is moving [13 . In other words, for very heavy 
quarks which slowly decelerate, the velocity is the only 
aspect of their initial conditions which infiuences the en- 
ergy loss rate. This turns out not to be the case for 
light quarks. Initial conditions for a classical string in- 
volve two free functions: the initial string profile and 
its time derivative. As we discuss in detail below, the 
instantaneous energy loss rate of a light quark depends 
strongly, in general, on the precise choice of these ini- 
tial functions. In the dual field theory, this reflects the 
fact that any complete specification of an initial state 
containing an energetic quark must also involve a char- 
acterization of the gauge field configuration. In the per- 
turbative regime, one can easily see that the interactions 
of heavy particles with a gauge field are spin-independent 
(up to 1/M corrections), but interactions of relativistic 
particles are spin-dependent even at leading order. So 
it is perhaps unsurprising that the energy loss of a light 
projectile also depends on the configuration of the glu- 
onic cloud surrounding the projectile in a non-universal 
fashion. 

One quantity which is rather insensitive to the pre- 
cise initial conditions of the string is the maximum dis- 
tance Axmax(^) which a quark with initial energy E can 
travel. It should be emphasized that we are consider- 
ing effectively on-shell quarks which can travel a large 
distance Ax before thermalizing. The maximum pene- 
tration depth Axmax grows without bound as the energy 
E increases. 

We numerically compute the penetration depth Ax for 
many different sets of string initial conditions, and find 
that the maximum penetration depth does indeed scale 
like E'^^^. Our results are illustrated in Fig. [l], where 
the logarithm of the penetration depth is plotted as a 
function of the logarithm of the initial quark energy for 
many different sets of initial conditions. As is evident 
from the figure, the penetration depth of a light quark is 
bounded by a curve Axmax = const, x E^^^. 

We also demonstrate the scaling relation Ax 

max ^ 

E^^^ analytically. As discussed in Ref. [34 , strings which 
correspond to long-lived massless quarks are approxi- 
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FIG. 1: A log- log plot of the quark stopping distance 
Ax as a function of total quark energy E for many falling 
strings with initial conditions of the form shown in Eq. (4.53). 



All data points fall below the red line given by Ax — 
{0.526/T){E/TVxy^\ 



mately null strings. A strictly null string is one whose 
worldsheet metric is everywhere degenerate. The qual- 
itative origin of this connection is easy to understand. 
Strings which correspond to light quarks fall into the 
event horizon. As they fall they become more and more 
light-like and hence closer and closer to a null configura- 
tion as time progresses. The profile of the null string is 
almost independent of the initial conditions used to cre- 
ate the string — for the quasiparticle excitations studied 
in this paper, the corresponding null strings are speci- 
fied by two numbers only, an initial inclination and ra- 
dial depth. By analyzing strings corresponding to light 
quarks as small perturbations away from null string con- 
figurations, we show that the total distance a quark 
can travel must be bounded by a maximum distance 
Axmax = {C/T){E/TVxy/^ for some 0(1) constant C. 
We numerically confirm that strings corresponding to 
long-lived light quarks are, in fact, close to being null, 
and obtain an estimate of the constant C. 

Although the endpoint motion of our string solutions 
is well approximated by appropriate light-like geodesies, 
consistent with the discussion of Ref. [34 , we find that 
the relationship between the parameters of the geodesic 
and the string profile and energy is more complicated 
(and rather different) than the surmises presented in 
Ref. [35]. This will be discussed further in Section [V|. 

In addition to studying the penetration depth, we also 
examine the instantaneous rate of energy loss, dE/dt. 
For light quarks the energy loss rate shows non-universal 
features and is sensitive to initial conditions. For the 
states we study, we find that it typically increases with 
time during the period when the dressed quark is a well- 
defined quasiparticle and sharply peaks during the final 
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thermalization phase. In other words, the thermahza- 
tion of hght quarks in strongly coupled SYM ends with 
an explosive burst of energy. This late time behavior is 
universal and independent of initial conditions. The fact 
that the light quark energy loss rate can increase with 
time is qualitatively different from the behavior of heavy 
particles [13 , whose energy loss rate monotonically de- 
creases. 

An outline of our paper is as follows. We define our 
conventions in Section [ill. In Section [Till, we discuss some 



of the subtleties involved in defining the light quark en- 
ergy loss rate and penetration depth, and spell out the 
relevant identifications between 5d gravitational and M 
field theory quantities. In Section [IV|, we discuss the 



dynamics of strings both from analytical and numerical 
perspectives. Implications of our results, and connec- 
tions with other related work, are discussed in Section [V|, 
which is followed by a brief conclusion. 



II. CONVENTIONS 




FIG. 2: Examples of perturbative diagrams contributing to 
the energy loss of a quark. One may regard time as running to 
the right. An energetic quark can scatter, emit gluons (which 
themselves may radiate or split into qq pairs), or annihilate 
with an antiquark in the medium. However, the total baryon 
number of the state remains constant. 



Five dimensional AdS coordinates will be denoted by 
Xm, while four dimensional Minkowski coordinates are 
denoted by x^. Worldsheet coordinates will be denoted 
as cr^ with a = 0, 1. The timelike world sheet coordinate 
is r = cr^, while the spatial coordinate is <j = a^. When 
discussing the dynamics of a single string endpoint, we 
will use cr* to denote the value of a at the endpoint. 

We choose coordinates such that the metric of the AdS- 
Schwarzschild (AdS-BH) geometry is 



ds' = 



-f{u) dt^ + dx^ 



du^ 



f{u) 



(2.1) 



where f{u) = 1 — {u/uhY L is the AdS curvature 
radius. The coordinate u is an inverse radial coordinate; 
the boundary of the AdS-BH spacetime is at = and 
the event horizon is located dX u = u^^ with T = (tti^/j,)"^ 
the temperature of the equilibrium SYM plasma. 



III. LIGHT QUARKS AND GAUGE/GRAVITY 
DUALITY 



Energetic quarks moving through a plasma are quasi- 
particles — they have a finite lifetime which can be long 
compared to the inverse of their energy. Some care is 
needed in defining the light quark penetration depth and 
the instantaneous energy loss rate. Fig. [2] shows some 
typical perturbative diagrams contributing to the energy 
loss rate of a quark. A energetic quark, scattering off 
excitations in the medium, can emit gluons which may 
subsequently split into further gluons or quark-antiquark 
pairs. The energetic quark may also annihilate with an 
antiquark in the medium. A natural question to consider 
when looking at Fig. [2] is which quark should one follow 
when computing the penetration depth? Once a quark 



has emitted a qq pair, or annihilated with an antiquark, 
it becomes ambiguous which quark was the original one. 
This issue is cleanly avoided if one focuses attention not 
on some (ill-defined) "bare quark", but rather on the 
baryon density of the entire dressed excitation. 

In QCD, or A/' = 4 SYM coupled to a fundamental rep- 
resentation M = 2 hypermultiplet, there is a conserved 
current which we will call ^baryon- Even though qq pairs 
can be produced by an energetic quark traversing the 
plasma, conservation of ^baryon ™plies that the total 
baryon number of the excitation will remain constant. 
The baryon density of an energetic excitation can remain 
highly localized for a long period of time. It is the col- 
lective excitation with localized baryon density which we 
will refer to as a dressed quark, or for the sake of brevity, 
simply as a quark. 

To evaluate the penetration depth of a quark, we will 
use the centroid of the baryon density. 



x{t) 



J d^x X p{t^ x) 
J d^x p{t^ x) 



(3.1) 



with p = ^baryon- "^^^ ccutroid x{t) givcs a natural mea- 
sure of where the quark is located at time t. At late 
times, when the quark has lost nearly all its energy and 
becomes thermalized, the dynamics of the baryon den- 
sity will be governed by hydrodynamics. In particular, 
at late times the baryon density must satisfy the diffusion 
equation 



{do-DV^)p = 0, 



(3.2) 



where D is the baryon number diffusion constant. When 
the diffusion equation is applicable, it is easy to see that 
the centroid ceases to move, dx/dt = 0. To define the 
penetration depth, we imagine measuring x{t) at some 
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early time t*. We then define the penetration depth Ax 
in the obvious manner as 

Ax = \x{oo)-x{U)\. (3.3) 

On the gravitational side of the gauge/gravity corre- 
spondence, the addition of a A/" = 2 hypermultiplet to the 
A/" = 4 SYM theory is accomplished by adding a D7 brane 
to the lOd geometry [38 . The D7 brane fills a volume 
of the AdS-BH geometry which extends from the bound- 
ary dit u = down to maximal radial coordinate 
and wraps an of the S^. The bare mass M of the 
hypermultiplet is proportional to 1/um [E], so for mass- 
less quarks the D7 brane fills all of the five-dimensional 
AdS-BH geometry. Open strings which end on the D7 
brane represent dressed qq pairs in the field theory. In 
the 5d geometry these strings can fall unimpeded toward 
the event horizon until their endpoints reach the radial 
coordinate Um where the D7 brane ends.^ For sufficiently 
light or massless quarks, Um > and open string end- 
points can fall into the horizon.^ 

The endpoints of strings are charged under a U{1) 
gauge field Am which resides on the D7 brane. The 
boundary of the bd geometry, which is where the field 
theory lives, behaves as an ideal electromagnetic con- 
ductor [39| and hence the presence of the string end- 
points, which source the D7 gauge field Am^ induce an 
image current density ^baryon boundary. This is 

illustrated schematically by the cartoon in Fig. [3]. Via 
the standard gauge /gravity dictionary [5l [H [71 [H [38] , 
the induced current density corresponding to each string 
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FIG. 3: A cartoon of the bulk-to-boundary problem at fi- 
nite temperature. The endpoints of strings are charged under 
a gauge field Am which lives on the D7 brane which 

fills the AdS-BH geometry. The boundary of the geometry, 
located at radial coordinate = 0, behaves like a perfect con- 
ductor. Consequently, the string endpoints induce a mirror 
current density on the boundary. Via gauge/gravity dual- 
ity, the induced mirror current density has the interpretation 
of minus the baryon current density of a quark. Similarly, the 
presence of the string induces a perturbation Kmn in the met- 
ric of the bulk geometry. The behavior of the metric pertur- 
bation near the boundary encodes the information contained 
in the perturbation to the SYM stress-energy tensor caused 
by the presence of the jet. 



^ One should bear in mind that even when the radial position of 
the string endpoints lies closer to the boundary than Wm, the 
string endpoints are nevertheless attached to the Dl brane, al- 
beit in the full lOd space. The embedding of the Dl brane is 
determined dynamically by minimizing the Dl worldvolume. In 
general, this means that the Dl brane wraps a 3-sphere inside 
the of the AdS-BH x background geometry. This 3-sphere 
varies in a non-trivial way as a function of the radial coordinate 
of the AdS-BH geometry. For a hypermultiplet with non-zero 
mass, the string endpoints must move on the internal as they 
fall down in the AdS-BH geometry, so that the string endpoints 
remain on the Dl brane. But for massless hypermultiplets, the 
corresponding Dl embedding is a simple product space, AdS-BH 
y^S^ . In this case, it is consistent to have the entire string sit at 
a fixed point on the while it falls in the AdS-BH background. 
Any additional motion of the string in the internal space will 
only add to the energy of the string without affecting its stop- 
ping distance and so will be of no interest for us — we want to 
find strings which carry a minimal amount of energy for a given 
stopping distance. In the large Nc limit, one can ignore the back- 
reaction of the D7 brane on the background geometry and the 
backreaction of the string on the D7, as well as potential insta- 
bilities involving string breaking or dissolving into the D-brane. 
These issues are discussed further in Section [V|. 

^ Strictly speaking, in the coordinate system we are using no por- 
tion of the string crosses the horizon at any finite value of time. 
Due to the gravitational redshift, the rate of fall du/dt decreases 
exponentially as one approaches the horizon. Nevertheless, it is 
natural to speak of the string endpoint falling "into" or "reach- 
ing" the horizon when u — Uh <^Uh- 



endpoint has a field theory interpretation as minus the 
baryon current density of a dressed quark. ^ 

The degree to which the baryon density is localized de- 
pends on how close the string endpoint is to the bound- 
ary of the bd geometry. The farther the endpoint is away 
from the boundary, the more the field lines of Am can 
spread out, and hence the more delocalized will be the 
induced image current ^baryon- limit where the 

radial coordinate U of the string endpoint is far from the 
horizon, U <^ u^^ the baryon density will be localized 
with a length scale ~ U [34 . We note that the appear- 
ance of the length scale U in the baryon density is natural 
since, for U <^ Uh^ it takes light an amount of time ~ U 



^ The fact that the induced mirror current density is minus the 
physical baryon current density is easy to understand. The 
baryon current density is given by the variation of the on-shell 
electromagnetic action with respect to the boundary value of the 
gauge field Am- The on-shell 5d electromagnetic action evalu- 
ates to a Ad surface integral, evaluated at the boundary with an 
outward pointing normal. In contrast, the image current den- 
sity induced on the boundary can be obtained by integrating 
the 5d Maxwell equations over a Gaussian pillbox which encloses 
the boundary. The resulting surface integral measuring the fiux 
involves an inward pointing normal {i.e., into the 5d bulk). 
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to reach the boundary. 

If at time the string's endpoint is at radial coordi- 
nate <^ then x{t^) approximately coincides with 
the spatial position of the string endpoint iCstring(^*) [34 . 
The string endpoint can only travel a finite distance be- 
fore falling into the black hole. The final spatial coordi- 
nate of the string endpoint iCstring(oo) will exactly coin- 
cide with ;^(oo) [34 . We therefore have 



and 



Ax |^^string('^) «^string(^*) I • 



(3.4) 



To make the quantity Axmax(£^) meaningful, we also 
need to measure the quark's energy at time t*. After 
all, we want to know how far a quark with a given ini- 
tial energy can travel. If the quark has been moving for 
some time prior to t*, it will have deposited energy into 
the plasma — we must disentangle the energy deposited 
in the plasma from the remaining energy of the quark 
itself. In the limit where the quark has an arbitrarily 
large energy which is localized in a arbitrarily small re- 
gion of space, separating the quark's energy from the 
energy transfered to the plasma will be unambiguous. 

Via Einstein's equations, the presence of the string will 
also perturb the bd geometry. As in the electromagnetic 
problem, the perturbation in the geometry will induce 
a corresponding perturbation in the Ad stress tensor on 
the boundary ^40^ ^ . The string itself has a conserved 
energy. For the states we consider in this paper, the 
endpoints of the string are very close to the boundary at 
time and hence have a very high gravitational potential 
energy. In Section |IV A| we argue that the energy con- 
tained near a string endpoint scales like l/u^. It is this 
UV sensitive part of the string energy that we identify 
with the energy of a quark. Via the gravitational bulk- 
to-boundary problem [also illustrated in the cartoon of 
Fig. ([3|] the high energy density near the string endpoint 
gets mapped onto a region of M space which coincides 
with the location of the quark's baryon density. There- 
fore, at time t*, we only need to compute the part of 
the string's energy which diverges in the ^ limit in 
order to identify the energy of the corresponding quark. 



IV. FALLING STRINGS 

The dynamics of a classical string are governed by the 
Nambu-Goto action 



Sng = —To / dr da 



-1: 



(4.1) 



where Tq = A/A/(27rL^) is the string tension, a and r are 
worldsheet coordinates, and 7 = det jab with jab the in- 
duced worldsheet metric. The string profile is determined 
by a set of embedding functions X^(r, a). In terms of 
these functions 



- 7 = (X • x'f - X^X 
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(4.3) 



lab = daX ■ dbX , 



(4.2) 



where X^ = drX^ and X'^ = d^X^ . 

The equations of motion for the embedding functions, 
as well as the requisite open string boundary conditions, 
follow from demanding vanishing variation of the Nambu- 
Goto action. Explicit forms of the resulting equations of 
motion, for the class of configurations we will consider, 
are shown below in Section IV A . The boundary condi- 
tions for the open string require that its endpoints move 
at the local speed of light and that their motion is trans- 
verse to the string. 

We will limit attention to configurations for which the 
string embedding only has non-zero components along a 
single Minkowski spatial direction which we will denote 
as X. We also restrict attention to initial conditions such 
that at worldsheet time r = 0, the string is mapped into 
a single point in spacetime. Explicitly, 

t{0,a)=tc, x(0,cr)=Xc, u{0,a)=Uc, (4.4) 

where the numbers tc, ^c, and Uc specify the 5d space- 
time location of the string creation event. The remaining 
initial data are the velocity profiles at r = 0, namely t, x 
and it as functions of a. One of these three velocity func- 
tions may be eliminated via gauge fixing. For example 
one may choose the gauge r = t. 

We will be interested in choosing initial data such that 
the subsequent evolution leads to configurations in which 
the two string endpoints propagate away from each other 
and become well separated before falling into the horizon. 
Choosing a frame in which the total spatial momentum 
of the string vanishes, this implies that one half of the 
string will carry a large positive spatial momentum in the 
X direction, while the other half carries a large negative 
spatial momentum. We also require that the velocity 
profiles are smooth near the string endpoints. A sufficient 
condition is that the Fourier series of the velocity profiles 
be rapidly convergent (pointwise). For brevity, we will 
refer to string configurations satisfying these conditions 
as "reasonable" . We postpone more detailed discussion of 
our specific choices of velocity profiles used for numerical 
studies to Section [IV Bl . 



A. Approximate solutions 

Consider a string created in the distant past. In par- 
ticular, take the radial coordinate of the creation event 
to be arbitrarily small, Uc 0. As time progresses, the 
string evolves from a point into an extended object and 
the endpoints of the string fall toward the horizon. An 
example of such string evolution (numerically computed) 
is shown in Fig. |4]. As we discuss in detail below, in the 
limit lie ^ 0, the string endpoints can be made to travel 
arbitrarily far in the spatial x direction before falling into 
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FIG. 4: A typical falling string studied in this paper, plotted in blue at four different instants in time. The string is created 
at a point and, as time passes, evolves into an increasingly extended object. Well after the creation event, but long before the 
plunge into the horizon, the string profile approaches a universal null string configuration which is largely insensitive to the 
initial conditions. Consequently, the string endpoint trajectories, shown in green and yellow, approach null geodesies. 



the black hole. Our strategy in this subsection is to con- 
struct an approximate solution to the string equations of 
motion which will provide a good description for times 
sufficiently long after the initial creation event but well 
before the string endpoints reach the horizon. This will 
be possible because, as we will discuss, at times well after 
the creation event but long before the final "plunge" , typ- 
ical string configurations approach near-universal forms 
which are characterized by only a few parameters. This 
observation will allow us to prepare states illustrating 
universal features and understand the resulting physics 
of quark energy loss, without requiring a detailed de- 
scription of the early-time dynamics responsible for the 
production of the quark- ant iquark pair. 

For reasonable falling string solutions, we will see that 
the endpoint motion is well- approximated by the trajec- 
tory of a light-like geodesic. Equations for null geodesies 
in the AdS-BH geometry are easy work out. For motion 
in the x-u plane, one finds 



dx 



geo 



dt 

du, 



geo 



dt 



f ie-f) 



(4.5a) 
(4.5b) 



where ^ is a constant which determines the initial incli- 
nation of the geodesic in the x-u plane and, more funda- 
mentally, specifies the conserved spatial momentum asso- 



ciated with the geodesic, f{u) dx^^^/dt = ^ 
over, we have 



More- 



We will be interested in string configurations where the 
spatial velocity of the string endpoint is close to the local 
speed of light for an arbitrarily long period of time (since 
this will maximize the penetration distance). Because 
open string endpoints must always travel at the speed of 
light, the velocity in the radial direction must be small 
and correspondingly, the radial coordinate of the string 
endpoints will be approximately constant for an arbitrar- 
ily long period of time. As the string endpoints become 
more and more widely separated, the string must stretch 
and expand. For reasonable string profiles, this implies 
that short wavelength perturbations in the initial struc- 
ture of the string will stretched to progressively longer 
wavelengths, resulting in a smooth string profile at late 
times. ^ Moreover, as the string endpoints separate, the 
middle of the string must fall toward the event horizon. 
This occurs on a time scale At of order u^. (This scale 
sets the infall time of a particle released at rest at the 
boundary, or of a null geodesic with ^ > 1.) 

The origin of this behavior can also be understood as 
follows. Consider the string at some time t shortly after 
the creation event. It will have expanded to a size ~ t. 
By construction, one half of the string will have a posi- 
tive large momentum in the spatial x direction, while the 
other half has negative x momentum. The spatial mo- 
mentum density must be highly inhomogeneous so that 
the two endpoints move off in opposite directions. As 
time progresses, the parts of the string with the highest 
momentum density will remain close to a string endpoint. 
Portions of the string with low spatial momentum den- 



dXcr 



du 



e-f 



(4.6) 



From this equation, one sees that geodesies which start 
close to the boundary, dit u = u^ ^0, can travel very far 
in the x direction provided ^ f{u^) ^ 1. In particu- 
lar, the total spatial distance such geodesies travel before 
falling into the horizon scales like u\/u^. 



"Unreasonable" string profiles can have structure on arbitrarily 
short wavelengths. While the initial structure will be inflated 
as time progresses, because the string endpoints can only travel 
a distance of order uf^/uc before reaching the horizon, one can 
always cook up initial conditions such that fluctuations in the 
string proflle never become small during this time interval. We 
will avoid such unreasonable initial conditions in this paper. 
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sity will lag behind the endpoints (in terms of motion 
in the x direction) and fall relatively unimpeded toward 
the horizon. Thereafter, the outer parts of the string will 
continue moving in the ±f direction while the endpoints 
slowly fall. This general behavior is clearly seen in Fig.|4]. 

With the above qualitative picture in mind, we now 
turn to the explicit analysis. To simplify the discussion, 
focus attention on one half of the string. We may choose 
worldsheet coordinates r = t and a = so that the 
embedding functions are determined by a single function 
x{t,u). The domain in the (t^u) plane in which x{t,u) 
is defined is bounded by a curve U{t) which defines the 
trajectory of the string endpoint. The location of this 
curve is fixed by the open string boundary conditions. 
With our choice of worldsheet coordinates, these bound- 
ary conditions are simply 



Gmn 
Gmn 



dX^ dX 



N 



dt dt 
dX^ dX^ 



dt du 



0, 
0, 



(4.7a) 
(4.7b) 



where the total time derivatives denote derivatives eval- 
uated along the curve U{t). As noted earlier, these con- 
ditions just express the constraints that the speed of the 
string endpoint equal the local speed of light with a ve- 
locity which is transverse to the string. 

With our choice of worldsheet coordinates, the deter- 
minant of the worldsheet metric is 



2/2 



[fx' 



/)• 



(4.8) 



Substituting this into the Nambu-Goto action (4.1), one 



finds the following equation of motion for the embedding 
function x(t, li). 



= 2u (1+/: 



/2\ 



2ti/ (/- 



Aufx'xx' 



+ 4/ [2 - / x' - 4{3-2f)x'x^ . (4.9) 

We want to construct an approximate solution for con- 
figurations where the string endpoint reaches the event 
horizon after traveling an arbitrarily large spatial dis- 
tance. We imagine first matching our approximate solu- 
tion onto an exact string solution at a time . At time , 
suppose that one endpoint of the string has fallen to the 
radial coordinate u^. Because we are considering Uc ^ 0^ 
we can always take 



Uc <^ <^ Uh . 



(4.10) 



For reasonable initial conditions, we can take suf- 
ficiently large so that the string will be close to a 
quasi-stationary configuration in which the string pro- 
file uniformly translates while the string endpoint slowly 
falls. In other words, we seek a perturbative solution to 



Eq. (4.9) of the form 



X{t, U) = Xsteady(^, u) + 5x{t, u) + 0{{5x) 



(4.11) 



where 



^steady (^, u) = (^t ^ Xq(u) 



(4.12) 



is a stationary solution to the equations of motion (4.9) 
and 5x{t, u) is a first order perturbation satisfying 



|<5i(t,u)| \5x{t,u)\<^\xo{u)\, 



(4.13) 



for all t>t^ and all u>U{t). 

At this point, the constant ^ appearing in the station- 
ary solution ^steady is logically independent from t he p a- 
rameter ^ characterizing null geodesies [c.f. Eq. (4.5)], 
but we will shortly see that the endpoint trajectory of 
the stationary solution ^steady is in fact directly related 
to the null geodesies discussed above. 

The endpoint trajectory may similarly be represented 
as a zeroth order curve plus a first order correction. 



U{t)=Uo{t)^SU{t). 



(4.14) 



where Uo{t) is the endpoint trajectory when Sx{t, u) = 0. 

The function 5x{t^ u) characterizes the perturbations in 
the string which have infiated to long wavelengths. Our 
basic strategy is to linearize the equations of motion and 
boundary conditions in both 5x(t^u) and 5U{t). From 
Eq. (4.9), the equation of motion for xq{u) is 







+ 4[(2- 



/)/-e'(3-2/)]^ 



(4.15) 



The general solution to this equation is given by functions 
which satisfy 



dxQ \ _ ^ 



(e - f) 



Cf) 



(4.16) 



where C is an integration constant. 

Neglecti ng th e perturbations 5x and 5U^ the boundary 
conditions (4.7) lead to the two endpoint equations 



dxQ 
du 



f 



fdu^y ie-f) 
\ dt ) e 



(4.17a) 
(4.17b) 



where all quanti ties are evaluated at th e string endpoint. 
Comparing Eq. (4.17a) with Eq. (4.16), we see that the 
two conditions agree provided C = 1. In other words, the 
boundary condition forces the integration constant C to 
equal unity. Furthermore, comparing Eq. (4.17b) with 



Eq. (4.5b), we see that the radial motion of the string 



endpoint in the stationary solution coincides with that 
of a light-like geodesic (when ^ of the stationary solution 
is identified with ^ of the geodesic). Since the speed of the 
string endpoint necessarily equals the local speed of light, 
this implies that the zeroth order endpoint trajectory, 
given by u = Uo{t) and x = A'o(t) = x steady {tMo{t)), is 
precisely a null geodesic. 
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With C = 1, the (negative root of the) differential 
equation (4.16) for xo{u) becomes 



dxo 
du 



f 



(4.18) 



(Taking the negative square root gives a solution which, 
for ^ > 0, trails the endpoint.) Su bstituting the steady 
state solution Xgteady into Eq. (4.8) and using the above 
differential equation for xo{u)^ reveals that the steady 
state string solution is one whose worldsheet metric is 
everywhere degenerate, 7 = 0. That is, Xgteady represents 
a null string which is everywhere expanding at the local 
speed of light. 

In the special case ^ = 1 (which will be of particular 
interest below), Eq. (4.18) may be integrated analytically. 
One finds 



xo{u) 



2 



tan 



^log 



(4.19) 



This is the well known trailing string profile of R ef. fT3^ . 
Similarly, when ^ = 1 the boundary condition (4.17b) 



may be integrated to find Uo{t). The solution is given 
implicitly by the equation 




FIG. 5: The inflation of a perturbation on a expanding 
string. The thin blue lines show the string at eight differ- 
ent instants of time. The uppermost black curve shows the 
endpoint trajectory. The perturbation to the stationary pro- 
file is the bump initially located close to the string endpoint. 
For clarity, we greatly exaggerate the size of the perturbation. 
The two red curves are lightlike geodesies which enclose the 
perturbation at all times. Even though the perturbation is 
initially highly localized, the two geodesies which bound the 
perturbation rapidly separate, and correspondingly the size of 
the perturbation rapidly inflates as it falls into the horizon. 



-xo{Uo) 



(4.20) 



But in much of what follows it will be useful to keep ^ 
arbitrary. 

Having found the zeroth order solution, we now turn to 
the first order correction which describes perturbations to 
a stationary null string. Linearizing Eq. (4.9) in dx{t^u) 
yields the equation of motion 



, 4C(C'-2/ + /2)^^ 2/2(1-2^2 + /) 



6x' . 
(4.21) 



A general solution to this equation can be constructed 
explicitly and has the form 

5x{t, u) = Uh Mz{t, u)) + g{u) ilj{z{t, u))] , (4.22) 

where (p{z) and iIj{z) are arbitrary functions, g{u) satis- 
fies 



and the function z{u^ t) is given by^ 



z(t^ u) 



[^steady (^, u) - Xgeo(^)] + Zq , 



(4.23) 



(4.24) 



^ The overa ll factors oiuu, n\lu^, and u^juj in Eqs. (4.22), (4.23) 
and ( |4.24| ) are inserted for dimensional consistency and later con- 
venience. 



with an arbitrary constant. For the special case of 
^ = 1, one has g(u) = —Uh/u. [Note that any constant 
of integration appearing in g{u) can be absorbed into the 
definition of ^{z).] 



Using Eqs. ( |4.6[ ) and (4.18), one easily finds that z(t^ u) 
is constant along null geodesies with the constant of mo- 
tion ^. Moreover, we may choose the constant zq such 
that z vanishes on the endpoint trajectory Uo{t). 

The fact that the perturbative solution to the string 
equations of motion contains two arbitrary functions (p{z) 
and iIj{z) is to be expected. As discussed in the previous 
section, the required initial data for the evolution of an 
initially pointlike string consists of two arbitrary velocity 
profiles. Evidently, the information contained in the ini- 
tial data gets mapped via the equations of motion onto 
the two functions (p{z) and tlj{z). 

It is easy to understand the physical nature of the per- 
turbations on top of the null string profile. The null 
string is everywhere expanding at the local speed of 
light. This expansion is analogous to cosmological infla- 
tion — perturbations defined on top of the null string 
at different points are causally disconnected and are 
transported along light-like geodesies. As illustrated in 
Fig. ([5|, neighboring geodesies increasingly deviate from 
each other. Therefore, as time progresses, the pertur- 
bations defined on top of the null string inflate to long 
wavelengths. 

To finish the first order analysis, we need to find the 
correction to the endpoint trajectory. At linear order 



in 5U and 5x, the boundary conditions (4.7) yield the 
constraints 



V-W = 0, 



(4.25) 
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and 
6U 



^(3/ -20 



(4.26) 

where (/?q = (^'{^)^ ipQ = '^'(0), and / and g are evaluated 
at UQ(t). 



1. Stopping distance 

If the size of the perturbations on top of the nuh string 
solution are small, then it is easy to compute the stopping 
distance — the total distance Ax traveled by the string 
endpoint after time t*. As the endpoint trajectories of 
null strings are lightlike geodesies, Ax is simply give n by 
the total spatial distance a geodesic travels. Eq. (4.6) 
gives the result for light like geodesies. 



dXg: 



du 



= ± 



(4.27) 



The ± sign reflects the fact that light-like geodesies can 
both fall toward the horizon and shoot upwards towards 
the boundary. Integrating this equation will yield the 
total stopping distance. 

Reality of Eq. (4.27) implies that > /(i^min), where 
Umin is the minimal radial coordinate achieved along the 
geodesic trajectory. We are interested in trajectories for 
which iXmin <C 1^*, so we require > /(iXmin) 1- Phys- 
ically, this corresponds to geodesies which only fall in the 
vicinity of and thereafter. This is sensible in the limit 
where the string creation point Uc ^ 0, since for reason- 
able initial conditions any portion of the motion in which 
the string endpoint is moving upward toward the bound- 
ary must occur during the initial transients shortly after 
the creation event and well before t^. 



For geodesies which only fall, Eq. (4.27) implies that 



Ax 



du 



ve-f{u) 



(4.28) 



With the restriction that ^ > 1, Ax is maximized at 
^ = 1. To lea ding order in u^/uh <C 1, evaluating the 
integral (4.28) gives 



(4.29) 



We must now relate the stopping distance to the quark 
energy at time t*. To make this meaningful, we want 
to estimate the minimum amount of energy required for 
a quark to travel a distance Ax before thermalizing. 
This requires determining how the string energy scales 
with u^. 

The canonical momentum densities of the string are 



given by 



TT 



M 




{X-X')X'^ -{X'^)X^ , (4.30a) 
(X • X')X^ - (X^) X'^] , (4.30b) 



where -7 = (X • X')^ - X^X'^. The energy of the string 
at time is then given by 



du \ 



(4.31) 



The zeroth-order approximation to the string solution 
is a null string, for which 7 vanishes. Hence, to describe 
a finite energy configuration, it is essential to include the 
perturbations to the null string profile. The determinant 
of the world-sheet metric will necessarily be proportional 
to the size of these perturbations. At linear order, and 
for ^ = 1, we have^ 



2L^ 



7 ^ 



(4.32) 



For a timelike worldsheet 7 must be negative, and hence 
we m ust have ?/;(z) < 0. Evaluating the string energy 
(4.31) to linear order in perturbations, one finds 



2tt 



j 

J U 



du 

U^f 



-2^lj{z{t,,u))]-^^^ . (4.33) 



This expression contains an infrared divergence near the 
horizon. This divergence refiects the unboundedly large 
amount of energy transfered to the plasma from the 
quark before time To extract a meaningful energy 
which can be associated with the quark at timet*, we 
focus on the UV sensitive part of the integral (4.33). Ne- 



glecting the IR region is tantamount to cutting off the 
radial integral at a radial coordinate iiiR < u^. The UV 
sensitive part of the string energy is the leftover part of 
the integral that diverges as -u* ^ 0. This is the portion 
of the string energy that should be identified with the 
energy of a localized quark jet at time t^. 

To minimize the energy (4.33) for a given value of u^^ 



one wants to maximize the magnitude of the function 
^l^{z) which characterizes the fiuctuation profile. How- 
ever, it is necessary to ensure that the perturbative treat- 
ment remains valid near the string endpoint — one can- 
not arbitrarily crank up the size of ^l^{z) as one must 



To linear order in the size of perturbations away from the null 
string, the function (p{z) appearing in Eq. ( |4.22| does not enter 
into the determinant of the worldsheet metric. In other words, 
perturbations in the steady state profile a^steady induced by (^{z) 
alone preserve (to first order) the everywhere null character of 
the string worldsheet. 

and 



More precisely, the upper limit of the integrals ( |4.3l| ) and ( |4.33| 
should not be w^, it should be the maximal radial coordinate of 
any point on the string at time — which rapidly approaches 
Uh- The contribution to the energy from the region w ^ 
refiects energy transferred to the plasma at times t <C • 
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ensure that the relations ( |4.13 ) are satisfied. More phys- 
icahy, we demand a weh behaved solution, which is ap- 
proximately a steady state solution, in the ^ 
limit. We remind the reader that, as discussed above in 
Section |IV A| , solutions which are approximately steady 
state solutions are dual to long lived quarks. 

At time t^^ the UV sensitive part of the string en- 
ergy comes from contributions near the string endpoint. 
Hence, we focus our attention on the region in which 
u — WU:^ with w — 0{1). Within this region we have 



w u% 



0{ul). 



(4.34) 



Then from (4.22) we see that 5x will scale with the same 
power of w« if 

</'W= (^)'<^W, (4.35a) 

^(^)= (j)'^(^), (4.35b) 

with (j){z) and il^{z) functions which remain bounded as 
0. With these scalings, a small 5x relative to xq can 
always be obtained by adjusting the overall normalization 
of 5x with a numerical factor which is independent of 
u^. Similar conclusions can also be reached regarding 
the scaling of 5U relative to that oiU.^ 

With the above scalin g of we see that the UV sen- 
sitive part of the energy ( 4.33| ) can be written 



ilVX 1 



where 
1 

C3 



TT 



^^^^ r 



-1/2 



(4.36) 



(4.37) 



and wiYi = uiYi/u^. By the scaling relations ( |4.35 ), the 
constant C is finite and independent of in the ^ 
limit. 



After using the result (4.36) to express in terms of 
E^, Eq. (|4.29D yields 



^^max(-^*) 



C f \ 



1/3 



(4.38) 



1 /3 

We reiterate that the EJ scaling is the maximum possi- 
ble power of energy consistent with the perturbative so- 
lution we have derived. In particular, it is the maximum 
power consistent with a string profile which is approxi- 
mately a steady state profile. 



From the differential equation ( |4.26[ ), the scahngs of Lp and -0 
imply SIA scales like w*. Therefore, in the ^ limit the 
smallness of 51A relative to lA, which is at most w*, can always be 
achieved by adjusting the size of 51A with a constant independent 
of w*. 



B. Numerical string solutions 

It is instructive to complement the above analytic anal- 
ysis with explicit examination of numerically computed 
string solutions. We wish to verify explicitly that {%) 
strings whose endpoints travel far in the Minkowski spa- 
tial are well approximated by null strings, (n) the end- 
point trajectories of such strings are well approximated 
by light-like geodesies with ^ = 1, and {iii) the maximum 
distance Ax that a string endpoint can travel scales like 
The first two points have already been demon- 
strated numerically in Ref. [34 , but we have extended 
that earlier analysis by exploring a larger sample of ini- 
tial conditions. This larger sample size is what allows us 
to address point {in). As we are using the same numeri- 
cal methods as in Ref. ^34j, this section closely parallels 
the analogous discussion there. 

To gain insight into the predicted E'^^^ scaling of Ax, 
we solve the string equations of motion numerically for 
a variety of initial conditions, and plot the penetration 
depth as a function of energy. As discussed below, we 



indeed find that the scaling relation (4.38) represents an 



upper bound on how far a string endpoint can travel for 
a given initial energy. Moreover, the numerical solutions 
provide a direct estimate of the constant C in the bound 
(|438l). 



For reasons discussed below (and earlier in Ref. p!3]). 
in our numerical analysis we have found it convenient to 
use the Polyakov string action. The Nambu-Goto action 
is classically equivalent to the Polyakov action 

Sp = -Y I d'ay^rj'^' daX^dhX^ Gmn , (4.39) 

where one has introduced additional degrees of freedom 
in r]ab^ the worldsheet metric. Varying the Polyakov ac- 
tion with respect to rjab generates the constraint equation 



lab 



\ rjab r]"""^ led ' 



This implies that 



-77 



(4.40) 



(4.41) 



so that the worldsheet metric differs from the induced 
metric only by a Weyl transformation. 



7?,5(r,a) = e2-(-'-) 7a6(r,a). 



(4.42) 



When Eq. (4.41) is substituted back into the Polyakov 



action, one recovers the Nambu-Goto action. 

The equations of motion for the embedding functions 
as well as the open string boundary conditions follow 
from variation of the Polyakov action with respect to the 
. Specifically, one finds 



da\ 



ab 



Gmn dbX^] 



ab 



dG 



NP 



dX^ 



(4.43) 
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together with the boundary conditions 
^Xf(r,a*) = 0. 



(4.44) 



Here a = a* denotes a string endpoint and tt^ is the 
canonical momentum flux on the worldsheet, 



8Sv 



N 



(4.45) 

We can fix the coordinate parametrization (r, cr) by 
choosing the worldsheet metric r]ab- As in Refs. [13^, ^S^l, 
we have found it convenient to choose rjab to be of the 
form 



Il^a5|| 



j:{x,u)-^ 



(4.46) 



We refer to S as the stretching function, which we take 
to be a function of x(r, a) and u{r^ a) only. The choice of 
the stretching function S is a choice of gauge. Changes in 
S lead to different embedding functions X^(r, a), but do 
not affect the geometry of the tar get w orldsheet. With 
a worldsheet metric of the form (4.46), the constraint 



equations Eq. (4.40) read 



X • X' = , 



(4.47a) 
(4.47b) 



Since we choose to study strings with point-like ini- 
tial conditions, the a derivatives are initially zero. 
Hence, we must choose initial time derivatives X^ which 



are consistent with th e constraint (4.47b) and the bound- 
ary co ndition (4.44). We may satisfy the constraint 



(4.47b) by fixing t in terms of x and u via 



ft 
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U 
J 



(4.48) 



To satisfy the open string boundary condition (4.44) at 
worldsheet time r = 0, we choose x and ii so that 



x'(0,cr*) =i/'(0,cr*) =0. 



(4.49) 



The set of pointlike initial conditions then reduce to the 



choice of two functions x and it obeying Eq. (4.49), to- 
gether with the initial radial coordinate Uc. 

To understand why it is preferable to start from the 
Polyakov action instead of the Nambu-Goto action when 
solving numerically for the string dynamics, note that 
the equations of motion (4.43) contain relative factors of 
{—r])~^ between different terms. Consequently, the string 
equations become singular whenever ^0. If we 

choose the worldsheet metric to be the induced metric, 
which is equivalent to starting from the Nambu-Goto ac- 
tion, then the equations of motion become singular as any 
part of the string approaches a lightlike configuration. 
This always happens at late times as the string acceler- 
ates toward the black brane. By using the Polyakov form 
of the string action, and exploiting the freedom to choose 



a worldsheet metric of the form (4.46), we may rescale 



the worldsheet metric so that the equations of motion 
remain well-behaved everywhere on the worldsheet. 

The energy of the string is a conserved quantity and 
can be computed from the data defining the initial con- 
ditions. With 



SSv 



St{r, cr) 



(4.50) 



denoting the conserved canonical energy density, the to- 
tal string energy is given by 



string 



d(j 7r^(0, <j) 



(4.51) 



Expressing this more explicitly in terms of the initial 
data, one finds that 



string 



27r i;(xc,iic)^c Jo 



dai{0,a). (4.52) 



1. Initial conditions and numerical results 

We consider a two-parameter family of initial condi- 
tions. Inspired by the strings studied in Ref. [34 , we 
choose 



x(0,cr) = A COS cr. 



(4.53a) 

^c\//(^c)(l-cos2a), (4.53b) 



and also take Xc = 0. As Uc ^ and Auc oo, these 
initial conditions generate strings whose endpoints travel 
arbitrarily far in the Minkowski spatial directions before 
falling into the black hole. Moreover, since ±(0, cr) is 
antisymmetric about cr = 7r/2 while u{0^ a) is symmetric, 
these strings are symmetric about x = at all times. 
These states therefore have zero total spatial momentum, 
but each half of the string has an energy and momentum 
that scale linearly with A for large A. 

As in Refs. [I3l|34], we choose a stretching function so 
that gradients of the embedding functions are small at 
all times during the string's evolution. We found by trial 
and error that stretching functions of the form 



S(x, u) 



1 



1 - u/uh 

1 - Uc/Uh 



<4.54, 



were adequate to generate long-lived strings with a vari- 
ety of initial conditions. In this work, the free parameter 
m was usually chosen to be 0.02. 



In terms of the initial conditions (4.53) and the stretch- 
ing function ( |4.54 ), the string's energy evaluates to 



E. 



string 



dG\/ A? cos^ cr - 



(1 -COs2cr)2 



27r 

(4.55) 

Because the string corresponds to a quark-antiquark pair 
in the dual field theory, £^string should be regarded as 
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FIG. 6: A plot of a numerically computed string at three different times, overlain with the analytic null string approximation 
(4.12). The string was created at a point at time t — with the initial conditions (4.53), for Uc — 0.014 and A — 2400. 
The corresponding energy E — ^string/2 85 700a/AT. The numerical string, shown as the solid blue curve, is plotted at 
successive times ti = 1-6/T, t2 — 11.4/T, and ts = 22.8/T, and the corresponding null string, shown as the dashed red curve, 
is plotted at the same times. The solid green and orange curves represent the numerically computed endpoint trajectories, 
and the overlain dashed blue curve shows the geodesic fit to the endpoint trajectory with geodesic parameter ^ = 1. The null 
string approximation agrees very well with the numeric string configuration at times t > a few Uh^ and the null geodesic curve 
likewise tracks the endpoint trajectory very accurately. 



twice the initial energy of a single quark. We empha- 
size that our strings are symmetric about x = O so that 
each half-string is approximately dual to a single dressed 
quark. In the following we therefore use 

E = \ ^string , (4.56) 

when discussing the dynamics of a single string endpoint. 

As in Refs. J3l [34], we used the Mathematica rou- 
tine N DSolve to integrate the string equations of motion 
(4.43) numerically. We chose the parameters A and 
in our initial conditions ( |4.53 ) from the intervals 

A e [1000,9000], (4.57a) 

— e [0.00867,0.019]. (4.57b) 

This generates half-strings with energies in the interval 



E 



e [26316,519031], 



(4.58) 



allowing us to study states with energies much larger 
than the characteristic scale a/AT (while simultaneously 
achieving high numerical precision). 

Our data runs stepped through this parameter space 
by fixing Uc and then generating strings for many val- 
ues of A. As a result, we were able to generate over a 
thousand string worldsheets and measure the associated 
stopping distances and energies. Our data are summa- 
rized in Fig. [l]. Each distinct line of data points in the 
plot comes from a single choice for Uc. 



2. Comparison to the approximate string solution 

Fig. |6] displays a typical numerically generated string 
at three different coordinate times. On top of the nu mer- 
ical string profiles, we also plot the null string (4.12) with 



^ = 1. Also shown in the figure are the endpoint trajec- 
tories overlain with the corresponding null geodesic with 
^ = 1. As is evident from the figure, the null string pro- 
vides an excellent approximation to the numerical string 
profiles for times t which are a few Uh or larger, and the 
difference between the actual endpoint trajectory and the 
null geodesic approximation is imperceptible. 

To further elucidate the quality of the geodesic approx- 
imation to the endpoint trajectory, we have computed the 
quantity 



E{t) ^ f(u(t)) {dx/dty 



(4.59) 



where X{ t) is the f -coordinate of the string endpoint. 
From Eq. (4.5a), one sees that for a geodesic S is constant 



and equal to ^. Over the course of the trajectory of the 
numerical string shown in Fig. [6], E! equals 1 to within one 
part in 10^, which is the limit of our numerical precision. 
Therefore, the endpoint path for this string is very well 
approximated by a ^ = 1 geodesic. We have verified 
similar results for many different sets of initial conditions 
which correspond to long lived quarks. 



3. Maximum penetration depth 

Our exploration of a wide range of initial conditions 
produced the data for penetration depths shown in 
Fig. [l]. All our data are consistent with the bound 



AXmax(^) 



0.526 



1/3 



(4.60) 



which explicitly confirms the scaling of the pene- 

tration depth derived in Section IV A . More generally. 



our numerical results clearly confirm the validity of the 
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asymptotic analysis leading to the approximate string so- 
lution presented in Section |IV A . 

It is instructive to estimate the maximum value of 
C b ased o n the perturbative analysis presented in Sec- 
tion IV A Clearly, from Eq. (4.37) one sees that C is 



maximized when the quantity 'il^{z) is maximized. How- 
ever, the validit y of th e null string approximation pre- 



IV A 



required <C 1. For 



sented in Section 

the integration appearing in Eq. (|4.37|) is order one. Fur- 



thermore, the value of C only depends on the cube root of 
the integral, so C is rather insensitive to its precise value. 
Therefore, in order to get a crude estimat e on the value 
of C we set the integral appearing in Eq. (4.37) equal to 
one. We therefore arrive as the estimate 



C 



V2 



: 0.45. 



(4.61) 



This is remarkably close to the numerically determined 
value of 0.526 

In addition to the sampling of point-like initial condi- 
tions yielding the data shown in Fig.[l], we have also stud- 
ied more complicated initial conditions describing strings 
which are not point-like at t = 0. The results obtained for 
these initial conditions also demonstrated the E^^^ scal- 



ing relation of Eq. (4.38), but generally yielded a slightly 



smaller value of C. All of our numerical results are con- 
sistent with the value for C determined from the data 
shown in Fig.[T], namely C = 0.526. However, we empha- 
size that this value, extracted from a finite sampling of 
initial conditions, is a lower bound on the true value of 
C. It is possible that a wider set of initial conditions will 
yield a larger value for C, although because of the close 
agreement with the estimate obtained in Eq. (4.61), we 



doubt that the true value is significantly greater than 
0.526. 



V. DISCUSSION 
A. Energy loss rate 

As Fig. [l] makes apparent, propagating light quarks 
in strongly-coupled A/" = 4 plasma do not have a unique 
stopping distance for a given energy. This result should 
not be surprising. Knowledge of the total energy (and 
momentum) of a quark-antiquark state is far from a com- 
plete specification of the initial state. The form of the 
disturbance in the gauge field (and other A/" = 4 SYM 
fields) will affect the subsequent dynamics. In the dual 
description, this additional information is encoded in the 
profile of the string connecting the quark and antiquark. 
Nevertheless, there is a rather simple characterization of 
the maximum penetration distance of a quark, scaling 
with energy as E'^^^. 

An interesting quantity to consider is the instanta- 
neous energy loss rate of a light quark. From the Ax ~ 
E'^^^ scaling of the penetration depth, one might expect 
that for light quarks the rate of energy loss per distance 



traveled, dE/dx (which essentially coincides with dE/dt 
while the excitation is a good quasiparticle), would scale 
like This expectation turns out to be incorrect, as 

we now discuss. 

Let /drag(^) — dp^/dt denote the four momentum lost 
by the quark per unit time. The long distance hydro- 
dynamic perturbation in the SYM stress tensor Tj^y^^^ 
is determined by the hydrodynamic constituent relations 
together with the energy-momentum conservation rela- 
tion ^31j . 



hydro 



(5.1) 



with F^{t,x) = -/LgW^^'H^-^quark(t)) the forCC- 

density (acting on the plasma) and iCquark(^) the quark's 
trajectory. 

As long as the quark's baryon density is well-localized 
in space, the energy loss rate may be determined by com- 
puting the energy fiux through a sphere Sr of radius R 
which encloses (nearly) all of the quark's baryon density. 
It is this energy flux which enters in the force density of 
Eq. (5.1). As 1/T sets the length scale on which a hy- 
drodynamic description of the stress tensor perturbation 
becomes valid [31], it is natural to take R ~ 1/T. The 
precise value chosen for R is irrelevant — during times in 
which the quark is a well deflned quasiparticle, its baryon 
density is localized over a scale <C 1/T while the distance 
traveled by the localized baryon density distribution is 
> 1/T. 

Using the dual gravitational description, one may com- 
pute the energy flux through Sr by solving the grav- 
itational bulk-to-boundary problem. Speciflcally, one 
solves Einstein's equations for the perturbation in the 
5d geometry due to the presence of the string and then, 
by analyzing the near boundary behavior of the met- 
ric perturbation [311 [411 [42] , extracts the change in the 
SYM stress tensor and uses this result to evaluate the 
energy flux through Sr. This procedure was carried 
out for heavy quarks moving at constant velocity in 
Refs. [311 EZl [33]. Carrying out the corresponding 
analysis for our non-stationary light quark worldsheets 
is computationally demanding and will be left for future 
work. However, there is a simple way of extracting the 
energy loss rate from the string proflle itself. The energy 
of a string is a conserved quantity. The high energy den- 
sity near the string endpoint should, as discussed in Sec- 



tion III , be regarded as the energy of the quark. This en- 
ergy is transported down the string by an energy flux tt^ 
[cf. Eq. (4.30b)] towards the event horizon. This energy 
flux corresponds to energy transferred from the quark to 
the plasma. Via the holographic bulk-to-boundary map- 
ping, this conserved flux is mapped onto the energy flux 
through Sr in the dual fleld theory. Without solving the 
bulk-to-boundary problem explicitly, one does not know 
a-priori precisely how to relate the bulk position {u^ t) 
at which one evaluates the energy flux down the string 
to a corresponding time and radius R of an energy flux 
measurement in the fleld theory. However, as long as 
the quark energy loss rate is changing sufliciently slowly. 
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FIG. 7: The instantaneous energy loss rate, dE/dt, of a 
highly energetic quark, normalized by its initial energy Eq. 
Instead of decreasing with time, as might have been expected, 
the light quark energy loss rate actually increases. At times 
near the thermalization time, which for this particular exam- 
ple is ttherm 24/T, the instantaneous energy loss rate grows 

like dE/dt ~ 1/V^therm - t. 



strings. The energy flux diverges for a null string. The 
flnite flux of the complete solution is determined by the 
function tlj{z{t^u)) [deflned in Eq. (4.22)] which charac- 
terizes the perturbation 5x{t^u) on the null string. This 
function is not universal, and depends on the initial con- 
ditions used to create the string. 

However, the late time behavior of the instantaneous 
energy flux is universal. As is evident from Fig. [7], near 
the thermalization time the energy flux down the string 
dramatically increases. This may be understood from our 
approximate string solutions. The energy flux down the 
string scales like (—7)"^^^ where 7 is the determinant of 
the worldsheet metric. For strings which are small per- 
turbations of null strings, Eq. (4.32) shows 7 is propor- 



tional to the function ip[z(t^u)) characterizing the per- 
turbations. Near the thermalization time ttherm = '^^/'^* 
and at a radial coordinate u corresponding to a flxed 
distance ~ 1/T from the string endpoint, the function 
z{t,u) behaves like 



u) 



t-t 



therm 



(5.2) 



retardation effects in the gravitational bulk-to-boundary 
problem can be neglected and the energy flux through 
Sr should be well approximated by evaluating the en- 
ergy flux down the string at a spatial distance ~ R from 
the string endpoint. 

In Fig. [7|, we plot the energy flux flowing down the 
numerically generated string shown in Fig. [6], evaluated 
at a distance of 1.75/(7rT) from the string endpoint. For 
this particular string, the string endpoint approaches the 
event horizon at a time t ~ 24/T, which should be 
regarded as the thermalization time ttherm of the light 
quark. As the string endpoint approaches the event hori- 
zon, the baryon density induced on the boundary rapidly 
spreads out and diffuses [34 . In the gravitational de- 
scription, this is due to the strong gravitational redshift 
incurred on the electric fleld sourced by the string end- 
points as they approach the horizon.^ As is evident from 
Fig. [7|, the energy flux down the string does not decrease 
in a power-law fashion as a naive T;2/^ scaling of dE/dt 
would suggest, but rather increases monotonically until 
the thermalization time! 

We stress that the precise form of the energy flux 
down the string is sensitive to the initial conditions used 
to create the string. This is easy to understand from 
the approximate analytic string solutions discussed in 
Section |IV A| . These approximate solutions, which cor- 
respond to long-lived quarks, are perturbations of null 



and hence becomes very small as t ttherm- By the open 



^ More precisely, as the string endpoints approach the horizon, the 
strong gravitational field of the black hole pulls the electric field 
lines, which are sourced by the string endpoints, towards the 
horizon. This results in the spreading out of the electric field 
lines and hence a spreading out of the induced baryon density 
on the boundary. 



string boundary conditions (4.7) and (4.25), the function 



^p{z) must vanish at the string endpoint which, as dis- 
cussed in Section [IV A| corr esponds to z = 0. Finiteness 
of the string energy ( |4.33 ) requires that '^'(0) be non- 
zero. Consequently, near the endpoint one may approx- 
imate iIj{z) ^ '^^ (Q) z. Neglecting the 0{u^/uh) correc- 
tions in Eq. (5.2), one flnds 



therm 



(5.3) 



We have numerically conflrmed the above scaling in the 

data shown in Fig. [7|. 

The late-time behavior ( |5.3[ ) implies that after travel- 
ing substantial distances through the plasma, the ther- 
malization of light quarks ends with an "explosive" trans- 
fer of energy to the plasma. This behavior is qualitatively 
similar to the energy loss rate of a fast charged particle 
moving through ordinary matter, where the energy loss 
rate has a pronounced peak (known as a "Bragg peak") 
near the stopping point. This peak in the energy loss 
rate has its origin in the energy dependence of cross sec- 
tions, which increase with decreasing energy due to the 
conformal nature of Coulomb interactions. 



In the gravitational description, the scaling (5.3) be- 



comes valid when the string endpoint starts to fall toward 
the horizon (i.e., when dlA/dX ceases to be small com- 
pared to one). As shown in Fig.[6], this happens relatively 
abruptly, so we expect the creation of large amounts of 
gravitational radiation to propagate to the boundary and 
induce a large perturbation in the SYM stress tensor cor- 
responding to this flnal burst of energy. However, we 
emphasize that because the energy flux flowing down the 
string is changing rapidly at late times, retardation ef- 
fects in the gravitational bulk-to-boundary problem can- 
not be neglected, implying that the result (5.3) for the 



15 



energy flux down the string cannot be directly equated 
with the fleld theory energy flux through a sphere Sr. It 
would, of course, be interesting to compute directly the 
energy flux in the plasma produced by the light quark 
jet as it thermalizes. Evaluation of the required bulk-to- 
boundary problem is currently in progress. 

It is interesting to speculate on the implications of our 
results for heavy ion collisions. Hard partons produced 
in the early stages of heavy ion collisions can traverse 
the resulting flreball and deposit energy and momentum 
into the medium. If the partons are moving supersoni- 
cally, their hydrodynamic wake will contain a Mach cone 
whose propagation can influence the distribution of par- 
ticles associated with a jet. If the hard parton under 
consideration is a very massive quark with mass m, the 
results of Ref. [13 predict an energy loss rate of the form 
dE/dx = dp/dt = —jj.p where, for strongly coupled SYM, 
/i = I \^T'^/m. Therefore, the energy loss rate falls ex- 
ponentially with time — heavy quarks in strongly cou- 
pled SYM lose the bulk of their energy in the early por- 
tions of their trajectories. The resulting sound waves, 
whose amplitudes will be largest at early times, may tra- 
verse much of the flreball before freezeout occurs. Con- 
sequently, sound waves produced by heavy quarks may 
be quite sensitive to medium effects and may experience 
substantial attenuation before freezeout. 

At least for strongly coupled SYM, the situation for 
light quarks is qualitatively different. As we have demon- 
strated in Fig.[7|, light quarks lose the bulk of their energy 
in the latter stages of their trajectories. The resulting hy- 
drodynamic wake will therefore have less time to attenu- 
ate and diffuse than is the case for heavy quarks. More- 
over, because the light quark energy loss rate increases 
with time, we expect the amplitude of the corresponding 
wake to also increase with time. Because of this, we ex- 
pect the spectrum and distribution of particles produced 
by light quark jets to be qualitatively different from the 
behavior of heavy quark jets. 



B. Fluctuations 

Throughout our analysis, we have treated the string 
dynamics classically. This approximation is valid in the 
limit of large 't Hooft coupling A. More precisely, a clas- 
sical treatment is valid in the limit that A ^ oo with 
E/^/XT flnite and flxed. (Recall that the string energy 
automatically scales like VX.) However, one would also 
like to understand when the classical analysis can be 
trusted if A is large but fixed. To determine this, one 
should compute the size of quantum fluctuations around 
the classical string proflle and compare the size of the 
fluctuations to the classical result. Natural speciflc quan- 
tities to consider are the fluctuations Ap in the quark 
momentum p. These fluctuations are deflned by the vari- 
ances 



If Ap is not small compared to p, then the reliability of 
the classical calculation is questionable. 

Formally, the mean momentum p and the connected 
correlator deflning (Ap)^ are both O(V^). Consequently 
\Ap/p\ = 0(A~^/^) and vanishes as A ^ oo. However, 
for large but flxed A the energy (and time) dependence 
of \Ap/p\ can be important. To see this, consider the 
case of fluctuations in the momentum of a heavy quark. 
Mean square fluctuations in the longitudinal and trans- 
verse components of the quark's momentum grow with 
time t as [15l [20] 



and 



(Apl)' - VXj^^^tH, (5.5a) 
(Apt)' - VXj^^^TH, (5.5b) 



respectively, where 7 = l/\^l—v'^ and v is the heavy 
quark velocity. Therefore, with large but flxed A, the 
relative size of fluctuations, Ap/p, becomes arbitrarily 
large both at sufficiently late times and (for longitudinal 
fluctuations) in the ultra-relativistic limit. 

To estimate the size of quantum fluctuations in the 
light quark's momentum, we will use the above results 
for heavy quarks as a rough guide. This is not unrea- 
sonable as the trailing string proflle used to compute the 
above momentum fluctuations coincides, in the v ^ 1 
limit, with the (^ = 1) null string derived in Section IV A . 
However, care must be taken — quantum fluctuations on 
top of the null string, which is a degenerate solution to 
the classical equations of motion, diverge. This is imme- 
diately apparent in the formulas (5.5) which blow up as 
V ^ 1. To estimate the momentum fluctuations for light 
quarks using these results, we must be able to associate 
the heavy quark velocity v (which is always less than 1) 
with the s ize of the classical perturbations Sx [deflned in 
Eq. (4.11)] to the null string proflle. To do so, we simply 
note that the trailing string proflle for a heavy quark with 
velocity v coincides with the null string proflle (which is 
the V ^ 1 limit) up to 0{l—v'^) corrections. Therefore, 
for the purpose of a rough estimate, we identify 
with the size of classical perturbations on top of the null 
string. As discussed in Section IV A 1 , if at time the 
radial coordinate of the string is ix*, then the largest^^ 
the perturbations to the null string can be is 0{ul/uf^). 
With the identiflcation 



0{ul/ui), we have 



\ 15/4 
J 



(5.6) 



During the portion of the quark's trajectory when it is 
a well deflned quasi-particle its momentum, by construe- 



(5.4) 



Strictly speaking, in Section |IV A l| we argued that Sx can be 
no larger than u^/uf^ only in the vicinity of the string endpoint. 
However, via the inflationary behavior of perturbations deflned 
on top of the null string (as shown in Fig. |5|, when t — = 
0(1/T) the perturbation in the string proflle will be determined 
by the near endpoint perturbations at time t = t^. 
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tion, is large and scales as p ~ ^/\u\/u^^. We therefore 
arrive at the estimate 



P 



9/8 



Tt 



Ai/4 



(5.7) 



In the limit u^/uh <C 1 {i.e. the high energy limit) we 
see that the size of quantum fluctuations are small rela- 
tive to the classical prediction for the momentum. But 
for fixed values of and A, the ^/i growth of the result 
(5.7) suggests there will be a time when quantum fluc- 
tuations become large. However, for the light quarks 
discussed in this paper, the quarks only exist as well 
defined quasiparticles up until the thermalization time 
^therm- As discusscd in Section pV A , in terms of the 
thermalization ti me is simply ttherm = ^h/^*- Evaluated 
at this time, Eq. (5.7) implies that quantum fluctuations 
are suppressed by a relative factor of {u^/uh)^^^ / X^^^. 
This is reassuring. 



C. Pair creation, string fragmentation and finite A^c 

In addition to large A, we have also assumed from the 
outset that A^c has been sent to infinity. This limit is 
what justifies the neglect of quantum fluctuations in the 
background AdS-BH geometry. At finite A/'c, string frag- 
mentation, backreaction of the string on both the geom- 
etry and the D7 brane embedding, and backreaction of 
the brane on the geometry also have to be addressed. 

Energetically, our string is unstable to breaking into 
many tiny pieces. This corresponds to quark-antiquark 
pair creation in the dual field theory. The small string 
coupling gs ~ l/Nc suppresses the amplitude for a string 
to break. Correspondingly, the rate of the string decay 
process is suppressed by l/N^, and hence the time scale 
for string fragmentation is, at large Nc, parametrically 
larger than any of the time scales considered in this work. 

Alternatively, the process of string fragmentation can 
also be described from the point of view of the D7 brane 
worldvolume as the spreading of a narrow flux tube, the 
original fundamental string, into more and more widely 
dispersed flux on the brane. In terms of the underly- 
ing string theory, the quanta of the worldvolume gauge 
field are little pieces of open string, so a uniform flux on 
the worldvolume is the same as a coherent cloud of lit- 
tle string pieces. Thinking of the dynamical instability of 
our string as a result of breaking into many pieces, or due 
to spreading into dispersed flux on the brane, are just two 
different descriptions of one and the same process which 
is suppressed at large Nc. 

Large Nc is also what justifies the neglect of backre- 
action of the D7 brane on the background geometry, as 
well as the backreaction of the string on the D7 brane 
embedding and on the geometry. Note that, as far as 
large Nc counting is concerned, the gravitational action 
scales as N^^ the action for the brane embedding and the 
worldvolume gauge field scales as Nc^ and the Nambu- 
Goto action describing the worldsheet of the string is of 



order 1. (In addition the three actions scale with the 't 
Hooft coupling as 1, A and a/A, but for now it is sufficient 
to focus on the Nc counting.) The brane is very heavy 
compared to the string, but still has a small tension in 
Planck units. Consequently, it is consistent to embed 
the brane in a fixed background geometry and then con- 
sider a string ending on the brane, without computing 
the 0{1/Nc) suppressed deformation of the brane which 
will be induced by the string. 

The issue of backreaction becomes more subtle when 
one solves for the linearized response of the metric in re- 
sponse to the string in order to determine the boundary 
stress tensor. The order A^^ stress energy of the string 
generates an order 1/N^ correction to the metric (since 
the 5d gravitational constant scales as 1/N^). Conse- 
quently, when evaluating the variation of the on-shell 
gravitational action, the perturbation in the geometry 
due to the presence of the string produces an order one 
contribution to the expectation value of the stress tensor. 

In addition to the string itself, another potential source 
for the stress tensor is the gauge field living on the brane 
which is sourced by the string endpoint. The 0(1) charge 
from the string endpoint gives rise to an order 1/A^c 
gauge field on the brane [as the gauge coupling on the 
brane is 0{1/Nc)]. Combining this with the overall Nc 
of the brane action would appear to give another order 
one source in the bulk, and hence another order one con- 
tribution to the expectation value of the stress tensor out 
on the boundary. 

However, it is important to note that the leading Am 
dependent term in the brane stress tensor is quadratic 
in the worldvolume gauge field, so that the order 1/A^c 
gauge field on the brane only produces a stress-tensor 
contribution of order 1/A^c- This conclusion is altered 
if there is an order one background electric field on the 
worldvolume (as in the original dragging string solution 
of Ref [13 ). In this case there are contributions to the 
brane stress tensor linear in the gauge field sourced by 
the string, which consequently give rise to an order one 
contribution to the stress tensor. Similar terms also arise 
if one studies finite, non-zero mass quarks where an order 
one background embedding scalar is turned on. Both 
of these contributions would need to be included if one 
wanted to generalize the stress-energy wake calculations 
of Refs. [30l [3TJ |32j [33] to finite mass quarks. 



D. Weak versus strong coupling 

It is natural to ask how the E'^^^ scaling of the penetra- 
tion depth in the strongly coupled limit compares with 
the analogous result for weakly coupled plasmas. When 
the 't Hooft coupling A (at scales ranging from the tem- 
perature T to the projectile energy E) is small, the en- 
ergy loss of a high energy parton moving through the 
plasma is dominated by near-collinear bremsstrahlung 
processes. The rate for an energetic parton (with en- 
ergy £^ T) to radiate a gluon which carries away 
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an 0(1) fraction of its energy, while interacting with 
a typical gauge field fluctuation in the plasma, scales 
as X^TI^/eJT gl SI SSI Si Si], up to factors de- 
pending logarithmically on the energy, which we ignore 
throughout this discussion. Therefore, the average dis- 
tance an energetic parton travels between emission events 
is Axrad(^) ^ ^JE/T/{\^T). The square root depen- 
dence on energy is due to LPM suppression, which is a 
consequence of multiple scattering during the formation 
time of a radiated gluon. 

Imagine creating a very energetic quark in a localized 
wave-packet with mean momentum p, and then measur- 
ing, at some later time, the total energy or baryon num- 
ber contained in a co-moving sphere of size R ^ 1/T ^ 
l/\p\ surrounding the wave-packet. The opening angle 
in near-collinear bremsstrahlung emission is parametri- 
cally small, AO - ^/X{T/E)^^'^. Therefore, the direc- 
tion of the leading parton is almost unchanged by these 
bremsstrahlung emissions. Since the speeds of ultrarel- 
ativistic excitations differ negligibly from the speed of 
light, this implies that all the partons produced by a 
cascade of near-collinear emissions have almost identical 
velocities. Consequently, near-collinear bremsstrahlung 
emissions do not significantly degrade the energy, or 
baryon number, contained in the co- moving sphere. As 
far as gauge-invariant measurements of energy or momen- 
tum are concerned, the entire collection of near-collinear 
partons behaves like a single collective excitation whose 
energy and momentum is nearly constant. 

This effective "quasi-particle" picture remains valid 
until the typical energy of the partons produced by the 
cascade ceases to be large compared to T. The typ- 
ical penetration depth will equal the radiation length 
Axrad{E) summed over the number of levels of showering 
which are required to degrade the typical parton energy 
from E down to ~ T. Since every emission transfers an 
0(1) fraction of energy to the emitted gluon, and every 
produced parton continues to shower, the typical energy 
of partons produced by a cascade with k levels of show- 
ering will be of order E/c^ for some c ~ 2. Therefore, 
the number of showerings required to thermalize an ex- 
tremely energetic parton grows only logarithmically with 
energy, and the total penetration depth differs from the 
radiation length for the first emission only by an 0(1) 
factor. The net result is that the penetration depth 
Ax{E) in a weakly coupled non-Abelian plasma behaves 
as ^/E/T/{X'^T) times factors depending only logarith- 
mically on E/T. 

Presumably, there is a smooth interpolation from weak 
to strong coupling in AT = 4 SYM. At intermediate cou- 
plings, the maximum penetration depth may be pro- 
portional to E^^^\ with an exponent i^(A) which varys 
smoothly from 1/2 as A Oto 1/3 as A oo. Alter- 
natively, the correct form might be a sum of two distinct 
contributions, T Ax = A{X) {E/TY^^ + B{\) {E/Tf/^, 
with A{\) = 0(A-2) and B{X) = o(A) as A ^ 0, and 
B{X) = 0(A-i/^) and A{X) = o(A-i/^) as A ^ oo. Sub- 
leading corrections to the weak-coupling energy loss rate 



which are suppressed by powers of the 't Hooft coupling 
are not known, and would be challenging to calculate. 
And subleading strong-coupling corrections, suppressed 
by inverse powers of A, are also unknown. Consequently, 
there is no way, at present, to determine a preferred in- 
terpolating form. 



E. Relation to other work 

In Ref. [35 , where the E'^^^ scaling was first proposed, 
various guesses for the analog of C were given based on 
different assumptions. The authors of this work were 
interested in calculating the penetration depth of a gluon, 
whose dual description was conjectured to be a folded 
string. The relevant string configuration was assumed 
to be given by a portion of the stationary trailing string 
profile of Ref. [13 , with the string (at any instant of time) 
coming up from the horizon, reaching a sharp hairpin at 
some radial coordinate u^{t)^ and then retracing the same 
path back down to the horizon. 

The authors of Ref. ^] estimated the penetration 
depth of a gluon by assuming that the hairpin in the 
string falls into the horizon along a lightlike geodesic. 
Without solving the string equations of motion, the pa- 
rameters of the geodesic were estimated in terms of 
and V. By relating the these parameters to the string's 
energy, the authors of Ref. [35] argued that the maximum 
penetration depth should scale like 



Axggpr 



Cggpr f E^ \ 
T \2TVxJ 



(5.8) 



The constant Cggpr was estimated to be between 0.35 
and 0.41. 

While we have found that the endpoint trajectories 
of strings corresponding to long lived quarks do follow 
lightlike geodesies (to quite high accuracy), the relation- 
ship between the parameters of the geodesic and the en- 
ergy of string are rather different than that presented in 
Ref. [35] . In contrast to the treatment of Ref. [35] , where 
the string energy was assumed to be well-described by 
Eggpr ^ VX/{u^Vl—v^)^ the energy of the strings con- 
sidered in this paper scale as \^u'f^/ul and the strings 
themselves are approximately null — this latter fact com- 
pletely fixes the corresponding geodesic parameter ^ in 
terms of the initial string profile via the equations of mo- 
tion. 

Despite these differences, it may be of interest to 
compare the penetration depth of a gluon estimated in 
Ref. [35 with the result for a light quark found in this 
paper. In doing so, it is natural to replace E^/2 E^, 



in Eq. (5.8) when converting from a folded string model- 
ing a gluon to an open string describing a quark. With 
this change, one may simply compare Cggpr to our mea- 
sured value of C = 0.526. Our result is larger than the 
estimates of Ref. [35^ by 30-50%. 

The E^^^ scaling of the penetration depth has also ap- 
peared in Ref. [36], which discussed the dynamics of jet- 
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like configurations in the bulk gauge field dual to the R- 
current in strongly coupled SYM. The coefficient of the 
scaling relation was not calculated in this work. How- 
ever, the coefficient characterizing i?-current jets nec- 
essarily differs from our result in Eq. ( 4.6Q| ), since the 
gravitational interactions of the 5d gauge field dual to 
the boundary i?-current are independent of A (at leading 
order in the strong coupling limit). 

We conclude our discussion by summarizing the 
physics which distinguishes light quark energy loss from 
that of heavy quarks. A comprehensive numerical study 
of heavy quark evolution has been performed in Ref. . 
Let us compare and contrast the behavior of heavy and 
light quarks. The penetration distance in both cases is 
non-universal for the same reason: the quark's evolution 
depends upon the initial gauge field. After several units 
of inverse temperature, the dual string in either case be- 
comes well approximated by small fluctuations on top 
of an analytic solution. As long as the quark is ultra- 
relativistic (regardless of its mass), the appropriate an- 
alytic string solutions are null strings, and the energy 
fiux flowing down the string is entirely determined by 
the non-universal small fluctuations. However, when a 
heavy quark has lost a sufficient amount of energy, its 
dual string profile will be well approximated by the non- 
null V < c solutions obtained in Ref. [13 . Thereafter, 
the heavy quark energy loss rate will be insensitive to 
fluctuations away from the analytic string profile, and 
the energy loss rate will simply be proportional to the 
quark's momentum. In contrast, the light quark energy 
loss rate remains sensitively dependent on fluctuations 
during its entire trajectory. This is a consequence of 
the fact massless quarks are always ultra-relativistic, so 
their dual string profile is nearly null at all times. Con- 
sequently, the energy loss rate profile of a light quark 
remains sensitive to the initial conditions for an arbitrar- 



ily long period until thermalization. 



VI. CONCLUSIONS 

Using gauge/gravity duality, we have studied the pene- 
tration depth of an energetic light quark moving through 
a strongly coupled A/" = 4 SYM plasma. An analytic 
asymptotic analysis shows that, for quarks which travel 
long distances through the plasma, the worldsheet of the 
dual string description nearly coincides with that of null 
string. Both the analytic analysis, and explicit numerical 
computations, show that for a given quark energy the 
maximum penetration depth Axmax(^) scales as E'^^^. 
Based on numerical results from a wide sampling of ini- 
tial conditions, we find Axmax(£^) = (C/T) {E/T^/Xf^^ 
with C ^ 0.5. We also find that the instantaneous en- 
ergy loss rate of a light quark is not universal. How- 
ever, independent of initial conditions, we find that the 
energy loss rate grows rapidly as the thermalization 
time is approached. Consequently, the thermalization of 
light quarks in strongly coupled A/" = 4 super Yang-Mills 
ends with an "explosive" burst of energy transfer to the 
plasma. 
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